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RADIATION PATTERN CALCULATIONS FOR AN ATS-F 
48-PANEL DEPLOYABLE-REFLECTOR ANTENNA 


R. Miezis 

Computer Applications, Inc. 
R. F. Schmidt 

’ Advanced Development Division 
Goddard Space Flight Center 


ABSTRACT 

Far-field radiation patterns are computed for a 
deployable -reflector antenna composed of N' panels 
which are, ideally, sections of parabolic cylinder. The 
vector Kirchhoff diffraction method (physical- 
optics) is applied at the antenna surface, and integra- 
tion is accomplished by digital computer to obtain the 
backscattered radiation patterns. Surface normals are 
found for indexed parabolic cylinders and distorted 
versions of this design structure approximating a 
paraboloid of revolution. It is shown that one of the 
methods evolved for a thermally distorted reflector is 
adaptable to the case when the assumption of straight 
lines on the metallized reflector mesh between support- 
ing ribs is not justified. 


GLOSSARY OF NOTATION 


Symbol 


E, H 


'P 


E r, E *, E * 

3 

' cr, £ 

F 

P 

p<t> h 


R 

a 


L 

Pi (y) 
A, A -1 
i. J, k 


u, v 
E, F, G 


Meaning 

electric and magnetic field vectors, subscripted 
unit normal to a surface 

scalar wave without time dependence, solution to wave equation 
radial distance 

spherical components of an electric -field 

electric-field of directive source, or a generic symbol for a function 
radial and angular variables of a cylindrical coordinate net 
focal length of a paraboloid 
position vector to surface 

partial derivatives of p with respect to subscripted variable 
focal length of a parabolic cylinder 

focal length of a parabolic rib lying in a parabolic cylinder 

antenna panel half-angle 

tangent to an arc or surface 

vector representing fabric line 

polynomial in variable y for i th antenna rib 

classical rotation matrix and its inverse 

Cartesian base vectors 

parameters for orthogonal net on parabolic cylinder 
functions of u and v 



GLOSSARY OF NOTATION (continued) 


Symbol 


Meaning 


n 

0, <I> 

& , 4 > 

0*0 

€ 0 

z o 

X, y, z 

x u* y u > z u 

x v> y v * z v 
d s 

L I 

D, R 

’ c 

k 

p e 

S 


symbol for intersection 

spherical angles for feed and surface discussion 

spherical angles for observer 

permittivity of free space 

permeability of free space 

intrinsic impedance of free space 

Cartesian components of a vector 

partial derivatives of x (u, v) with respect to u, etc. 

partial derivatives of x (u, v) with respect to v, etc. 

differential area 

sampling interval in numerical integration 

diameter and minimum radius of curvature of scattering object 

wavelength 

feed position vector 

source strength of a feed element 
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RADIATION PATTERN CALCULATIONS FOR AN ATS-F 
48-PANEL DEPLOYABLE-REFLECTOR ANTENNA 


INTRODUCTION 

Consideration of continuous reflecting surfaces with discontinuous deriva- 
tives for the Applications Technology Satellite (ATS-F) has led to the develop- 
ment of methods for computing normals, areas, etc. of such surfaces as re- 
quired by the complex-vector Kirchhoff diffraction theory. Departure from the 
class of surfaces obtained by rotating arcs, derived by conic sections, has re- 
sulted in methods which are partially analytic and partially numeric in nature. 

For example, the normals to a parabolic -cylinder section can be obtained 
analytically, but the normals for 48 indexed sections of this type are injected 
into the program using the analysis for a single panel in conjunction with a 3 x 3 
rigid-body rotation matrix. Actually, the mixture of analytical and numerical 
techniques is also seen with surfaces obtained through conic section, although 
to a lesser degree* In any event, the reader will not find a single equation, or 
even a set of equations, which can be combined to form a "master equation" for 
these diffraction problems. Rather, a process, is developed, whereby the quan- 
tities required by the Kirchhoff formulation are made available to the computer 
for numerical integration. 

A brief review of the mechanics of the Kirchhoff method is given here; the 
referenced documentation provides greater detail. In the far-field (Fraunhofer 
region) the general formulation for the electric field reduces to a single integral 
which contains a mathematical simulation of a sheet current that flows tangentially 
along a reflector surface. 


E (x\ y', z') *-l£L!L 

4 rr 


f n x H 

Js, 


i^s 


It can be seen that "gain" is implicit in the parameter co (angular frequency). 
Also, it is noted that the vector cross product of the unit surface normal n with 
the magnetic field Hj necessarily results in a sheet current flowing along the 
local tangent to the surface. This is consistent with the physics of the situation 
and agrees with the intuitive concept of a current on a reflector antenna. Under 
the integral sign, the factor 0 = e' ikr /r is recognized as the product of a phase 
term, the periodicity of Young, and a space attentuation or divergence term. 
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Annexing differential area, ds, completes the simulation. The result is an inte- 
gral over a collection of Huygens sources, in vector notation, for three-dimen- 
sional bodies that are perfectly conducting. 

The preceding views the Kirchhoff formulation from the location of the 
observer of a field distribution. The calculation of the magnetic field distribu- 
tion H t at the reflector is an important preliminary detail which can be some- 
what involved, depending on the feed complexity. Feed directivity, polarization, 
space divergence and periodicity between feed and reflector, and reflector 
boundary conditions are taken into account to obtain the distribution from the 
incident magnetic field H. . The relationship between the incident field and the 
distribution H x which is integrated is 


mm ** 

n x Hj s 2 n x Hj 


The boundary condition, therefore, enhances tangential H i and annihilates normal 
H i . It is sufficient Jo calculate only the scattered fields (x'.y'.z') in theFraun- 
hofer region since H (x',y',z') differs from the former by only a physical constant 
z o “ Mq /£ 0 anc * wi orthogonal rotation in space. These scattered electric and 
magnetic fields are synchronous in time (cophased). 

In summary, the vector Kirchhoff formulation deals with surfaces, and 
integration is over the actual antenna surface rather than over an ’'aperture- 
plane.” The simulation closely adheres to the actual physical process involved 
in the illumination of the antenna by converting incident fields to electric sheet 
currents. Polarization, amplitude, and phase are carried throughout the entire 
problem, which is distinct from ray-optics analysis. Finally, since the present 
simulation deals with complex vectors, the- results can be regarded as three 
spatially orthogonal complex- scalar solutions for each problem, of which only 
the transverse fields and are retained for the far-field, together with 
their respective phases. 


ANALYTICAL APPROACH TO SURFACE NORMALS 

The unit normals to a given surface can often be obtained in several ways, 
however, as the description of the surface becomes more tenuous some of these 
methods cease to be applicable. Prior to dealing with highly arbitrary surfaces, 
it is helpful to review a few standard approaches for well-known geometries , v; 
The subsequent transition to more general cases then becomes obvious. In the 
following discussion, no attempt is made to obtain an analytic "best-fit" surface 
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in the BMS sense, or any related approach, in order to obtain the normals. In- 
stead, the goal is simply to provide the computer with the Cartesian components 
( n x> n y , n t ) of a unit normal for each curvilinear quadrilateral (differential 
area) where the surface may be specified by only a limited amount of data and 
perhaps some auxiliary constraint or assumption such as straight fabric lines. 

The paneled antenna considered for ATS-F is a close approximation of a 
paraboloid. The latter is described by 


xsc sin £, 


y 


<J COS £, Z 



(Zj a constant), 


and this parametric form will serve to introduce methods of computing normals. 
Another way of describing the same paraboloidal surface is 


- cl (x, y) = : 


x 2 + y 2 


TF 


+ or 3 2 (x, y, 2 ) =3 1 (x, y) - z- 0. 


In one text 1 , the unit normals to a surface are found using 


- / 3 z f Bz* t 

n s: ( “• y j* *• J + k 




+ 1 


i/2 


This method is a variation of the so-called gradient method. 2 The gradient 
method, in general, terms, leads to the form 


V3 0 (x, y, z) 


ns* 


V2L (x 


.y, z) / ,33, t 33 ,33 2 \/f/B3 2 \ 2 /^3 \ J /33 \ 2 "1 ,/2 

—\ t 1 t?-> of w/lw) Aw) * \w / J • 




/ 2 


1 Ref, 1, page 208 

2 Ref. 2, page 27 and Ref, 1, page 103 
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Another method of obtaining unit normals to surfaces is suggested by con- 
sideration of a pair of tangent lines at a point . 1 A normalized vector cross- 
product of the tangents is 


n = P<r * Per x \\ 


where p = I x + J y + kz is the vector from the origin to the surface, and the 
subscripted p denotes partials with respect to the indicated parameter. The 
reader can verify that all three methods lead to the same result for the family 
of paraboloids given above: 


n * (-Icr sin £ + J crcos £ + k2F)/ (cr 2 + 4F 2 ) 17 2 


It is noted that all of the preceding methods are general in the sense that 
they can be applied to other parametric representations 2 for surfaces. A de- 
parture from surfaces of revolution will now be made, with the final objective 
of extending the technique to distorted surfaces and continuous surfaces with 
discontinuous derivatives. 

The idealized paneled structure, approximately a true paraboloid on ATS-F, 
consists of indexed sections of a parabolic cylinder whose focal-length is taken 
as F (a line-focus). The normals to a parabolic cylinder described by 


« cr 2 

x - cr, y - any value, z = _ — + zi, (z, = constant) 

4 F. 1 


can be found using a variation of the gradient method since 


z = 3 i (x) = 4T + Z > 

C 


1 Ref, 1, page 209 

2 Ref. 3, page 13 


4 



now represents the surface. Then 


n = (- cr I + 2 F c £)/( a 2 + 4 F 2 c ) l/2 


for the family of parabolic cylinders. 

The crossed -tangent method can also be applied to the parabolic cylinder. 
For the case of the paraboloid two orthogonal 1 tangents, p„ and p r , were used. 
Here, two orthogonal tangents, p and p , serve the same purpose. 

x y 

P = i x + f y + £ + z\, (7 X = constant) 


Then 


and 




which are obviously orthogonal, give the same n for the parabolic cylinder via 


Py 

Attention is called to the fact that normals to surfaces are sought for the 
vector Kirchhoff formulation. The selection of simple curves (Jordan arcs) lying 
in a given surface allows the computation of a disc of normals at each point on a 


1 The tangents for the paraboloid are obtained from an orthogonal net (or, Q. Verification of 
their orthogonality is seen in the vanishing dot product. 

P<r ' Pt, - ' s ' n £ “ i cos £ + £ ^ * 0 V cos C + f cr sin Q * 0. 
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curve. It is usually not obvious which of this infinity of normals at a point on the 
curve is also the unique normal to the associated surface. The risk of making 
an error here is obvious since at least two distinct curves are required through 
a common point before a surface normal can be defined. The intersection of two 
discs of normals is the desired surface normal. For the preceding reasons, 
single arcs will not be used here to infer surface normals. 

. * 

The surfaces considered here are all simple and have unique normals with 
but one exception. A pair of normals exists at every point where panels are 
contiguous, corresponding to the normals of the "parent’' parabolic cylinders. 
Practically, this causes no difficulty with the Kirchhoff integration in the pro- 
gram since normals are always taken in the interior of curvilinear quadrilaterals 
(differential areas). Implicitly, however, this represents a departure from rigor 
since the Green's theorem underlying the Kirchhoff formulation assumes con- 
tinuous functions. In view of other approximations such as quantization of the 
illumination distribution, and the use of limited data on the surface, the de- 
parture is not considered serious. It is noted that a finite or mathematically 
open surface also represents a discontinuity. The Kottler contour integral 
which compensates for this is not considered here as the transverse far-fields 
are not affected. 


NUMERICAL APPROACH TO SURFACE NORMALS 

A somewhat more realistic simulation of the ATS-F deployable antenna is 
obtained by admitting distortion of the parabolic i ribs supporting the (metallized) 
reflecting mesh. If it is assumed that the reflecting mesh is defined by straight 
lines (L) between supporting ribs 2 , as shown in Figure 2, there are several ways 
of obtaining surface normals. All of the methods considered at this writing are 
predicated on the formation of a polynomial z = P(y) for each supporting rib. It 
is anticipated that the number of data points describing a rib condition will be 
much smaller than the number of samples required in the radial direction to ac- 
complish the integration of the Kirchhoff formulation. The discussion is re- 
stricted to a single panel throughout since the treatment for each panel is identi- 
cal and the indexing of panels, to obtain the proper resolutes of the surface 
normals in 3 -space, follows from simple rigid-body rotations of classical me- 
chanics. 3 A Cartesian (x, y) net replaces the cylindrical (cr,£ ) net used with 
axially symmetric surfaces. 

Appendix A. It is shown that the intersection of a plane and a parabolic cylinder (F ) is a 
parabolic arc (F R ) or rib, F R =F c /cos 2 a, where ais the panel half-angle. 

The straight-line assumption will be removed in the subsequent development. 

3 Ref. 4, page 107 
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z 



Figure 1. Paneled reflector with distorted ribs subdivided 
by means of a system of (x, y) coordinate nets 


When rib coordinate data is given in Cartesian (x, y, z) form, a vector with 
origin atO represents a point on the rib and, therefore, on the panel. Rib data 
can be transformed 1 to the (-y, z) half -plane via the indirect transformation 
(A* 1 ) so that z = P(y), a polynomial in y alone, can be formed. This is effectively 
a rotation of the rib about the z-axis. If this operation is performed for all the 
ribs a number of elementary techniques 2 , such as the direct-transformation (A), 
can be used in conjunction with the (x, y) net to provide rib coordinates (x, y, z) 
to as fine a resolution, or sampling, as may be required by the Kirchhoff method 
at a given frequency. Once the rib coordinate data (between given data points) is 
available, it is possible to write the components of the taut-mesh line L as the 
difference cf rib position vectors Vj andV 2 . See Figure 1. This line is one of 

1 Ref. 4, pp. 108-109 describes the tronsformotions A and A* 1 . 

^ The details of the process are contained in the program listing. 
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the tangents to the distorted surface under the straight-line fabric assumption, 
and immediately suggests the first method for finding surface normals via an 
interpolation process. 


An obvious means of obtaining a normal to the surface at a panel edge is to 
obtain another tangent somewhere, and use the method of crossed tangents, 
properly normalized. The polynomial Pj (y) for rib R, can be differentiated to 
obtain the slope, r x , everywhere along that rib. Similarly, P 2 (y) for rib R 2 
will supply another tangent, , but the latter vector must undergo a rotation 
transformation about the z-axis in the amount 2a, the panel angle. Tangents to 
the polynomials necessarily have a zero component in the x-direction since 
(y) is a function of y alone. It can be seen that the rib tangents are not, in 
general, orthogonal to the taut-mesh line L . This does not preclude their appli- 
cation in the crossed-tangent technique for finding the normals to the surface. 

It is sufficient that two distinct tangents lying in a tangent plane have been found , 
and these uniquely determine the normal in at a point on the surface. With a 
normal rij and n 2 at either end of a strip of fabric, as shown in Figure 1, other 
intermediate normals are found by simple linear interpolation of the difference 
of the components of and n 2 along L = | L | . As a special example, the dif- 
ference 


( n xl - n x2> 1 + ( n yl - n y2> 5 + ( n *l - n z 2 > * = < 0 ’ °’ 0 > 


for the case of the parabolic cylinder panel so that no incremental correction 
to the normal at n 1 is applied while traversing a strip of length L. The normal 
is everywhere the same between points 1 and 2, as it should be, for the indexed, 
ideal parabolic cylinder. 

A second method of obtaining the normals to the distorted surface can easily 
be evolved from the preceding. The taut-fabric concept that regards the curve 
between points 1 and 2 as a perfectly straight line, mathematically, can now be 
removed. Since two tangents suffice to form the normal, three closely-spaced 
points lying on the surface of interest lead to a normal. If the surface is ade- 
quately specified and subdivided, by an (x, y) net here, three of the four coordinate 
triplets defining an increment, of area can be used to form a pair of difference 
vectors (i.e., tangents). Interpolation of the normals is not necessary by this 
method. The surface may be specified by a large amount of data on the rib and 
fabric positions or, alternatively, by a description of the curves between ribs 
(points 1 and 2). The latter may be empirically derived. Should the fabric hap- 
pen to form straight lines as before, an interpolation, not of the normals to the 
surface, but one based on the z-coordinates at the end points of a strip can be 
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made. The 3 -point, two-tangent approach for obtaining normals is obviously not 
restricted to an (x, y) net or panels, etc., but is generally applicable. 

In conclusion, analytical methods are difficult to apply where surfaces are 
specified by data and auxiliary rules or equations. Methods which utilize the 
given data directly, or introduce the intermediate step of forming polynomials 
to obtain additional coordinates between data points, are therefore employed in 
the absence of an analytical surface description. A comparison of three distinct 
subroutines, representing varying degrees of departure from a parabolic cylinder, 
has been presented here: 

(1) ideal parabolic-cylinder panel 

(2) linear-mesh panel with distorted parabolic ribs 

(3) arbitrary panel surface with distorted ribs. 

An assumption throughout has been that the antenna has piece-wise (panel-wise) 
continuous derivatives. 


INCREMENTAL AREA 

The computation of incremental area for the paneled antenna is closely bound 
to the subject of surface normals. When analytical surface descriptions are 
available standard 1 methods apply. For example, the ideal parabolic cylinder 
panels can be treated by letting 


so that 



E = x u + yi + z u ~ 1 + — » F = 0 (orthogonal net), G = x 2 + y 2 + z 2 = 1 

4 F 2 v v v 

C 

and 2 


d S = (E G - F 2 ) 1/2 du dv = (y 2 - yi ) (l + JlMdx. 

\ 4 F 2 / 

C 


1 Ref. 3, p. 13 

2 Ref. 1, p. 206 
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Then an incremental area can be evaluated as 



d s = s =. (y 2 - y]l ) 


x 



Vx 2 + 4 F 2 +F In 

c c e 



7x 2 + 4 F 2 

c 

2F 

C 



Several general expressions 1 can be written for ds, independent of any 
specific definition of surface. Using 


" = (£ u x P v )/( EG-F 2 )*/ 2 , 


d s _(EG- F 2 ) 1/2 dudv = I P u x p v | du dv = sec /n, k du d v. 


All three can be applied to the ideal parabolic cylinder, and lead to 


sec 




The last method is also applicable when distorted surfaces are encountered 
since the increment of area on the curved surface is related to the increment 
of area on the chosen planar net, the (x,y) net here, by a simple projection. 
Normals to the distorted surfaces are found by methods such as the non- 
orthogonal crossed-tangent approach, described previously. Then the angle 
between the base vector k and the normal to each increment of area can be 
computed and used to evaluate ds. 


As » sec /n, k A A 

(curved surface) (planar-net area) 


1 Ref. 1, pp, 208-209, Ref. 5 p. 48, Ref. 6, pp. 321-331, Ref. 7, pp. 81-83 
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The computation of A A is simple and depends on the net selected for the problem: 
(x,y), (cr t 0, etc. 


NUMERICAL INTEGRATION 

The preceding sections dealt with the normals and the differential areas 
required by the Kirchhoff formulation, and these are now incorporated in die 
numerical integration. Since the original form for the backscattered field, 


E (x', y', z') - - 1 f n x H, 0 d s 

4 n Js, 

cannot be evaluated directly by purely analytical methods, in general, it is 
replaced by 


§(x', y', z') 



The continuous approach is, therefore, replaced by a discrete or quantized ap- 
proach. For this reason no "master equation" is set down for the solution. The 
quantities , and AS^ are properties of the specified surface: normals and 
curvilinear quadrilaterals of area. Illumination distribution and the phase- 
factor 0^ are determined by selecting the radiation source and the locus of 
pattern observation. 

In the present ATS-F problem, a summation of the type shown above was 
performed for each of the 48 reflector panels. An (x,y) orthogonal net was con- 
venient for subdividing the surface. There is no restriction on the method, how- 
ever, and other geometries can be treated by the same process. Whenever pos- 
sible, a coordinate net that is natural for the given geometry should be employed. 
Paraboloids for example, are most easily treated by a cylindrical or (cr, £) net, 
but this would be a poor choice for the ATS-F paneled reflector problem. Since 
the latter exhibits straight panel edges, an error in total area would be introduced 
by the (cr, £) net unless special devices were employed to delete certain (cr, £) 
quadrilaterals. At best, a residual error due to a saw-tooth edge would result. 
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The question of adequate sampling 1 of the surface arises since the continu- 
ous and discrete forms for E (x, y! z ) are equal only in the limit as As ^ becomes 
arbitrarily small. In the computer program an integration control constant 
designated LI (any positive real number) is made available to the program user. 
This number defines the maximum length of an edge of a curvilinear quadri- 
lateral in wavelengths. For the present problem A x « Ay = LI k , and LI = 0.5 
or 1.0 is usually adequate. The stability of solutions with respect to the sampling 
interval can be verified occasionally by computing a few field points at LI - 0.25, 
particularly for large observation angles (9) from the reflector-system axis. 

A program-related topic of great practical significance is the fitting of poly- 
nomials to given rib data . The least squares method was used to establish the 
coefficients of z = P(y) up to and including y s . In order to verify the fit, an .RMS 
figure is published by the program for each rib, and an overall RMS figure is 
provided for the set of ribs. For the present problem order 2 yielded the best 
fit for thermally distorted rib data at HOUR 11 and HOUR 4, the RMS value being 
less than *002 inch for the set of 48 ribs. The worst case observed yielded .007 
inch RMS for one rib; the exact reason for the large RMS is unknown. Most in- 
dividual ribs exhibited a polynomial fit well below the desired .005” criterion 
established by the ATS-F Project Manager for the 30-foot reflector. 


PROGRAM INPUT-OUTPUT QUANTITIES 

Program input can be separated into four groups as follows. 

(1) Option Control Constants 

These constants specify the number of feeds, the type of reflecting surface 
(paraboloid, paneled reflector, etc.), and the number of ribs, where applicable. 

(2) Reflector Inputs 

The parameters entered here are: maximum aperture radius, minimum 
aperture radius when a central region of the reflector is to be deleted, distance 
of the reflector "vertex" from the origin of coordinates, and focal length when a 
paraboloid or parabolic -cylinder paneled structure is selected. Data may be 
supplied in lieu of an analytical surface description. 


1 Ref. 8, pp. 321*323. The author discusses an antenna analogue to the classical sampling 
theorem of information theory by Shannon, 
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(3) Feed Inputs 


Fox’ each of the specified feeds, the source strength, initial or starting phase, 
polarization vector components, directivity, and feed position/orientation are 
read in. The associated symbols are Sj,^ , 5 . orh l , N in cos N 0 , , and 

Euler angles a, f , y. The operating frequency is also included here. 

(4) Locus of Observation 

Several choices can be made for the locus of observation of the radiation 
pattern. The usual one is a spherical locus given by range r = constant, angle 
4> = constant, and an angle $ « variable. Here $ <, $ < n , and the in- 

crement A 0 determines the number of coordinate points at which the field is 
computed. 

Program output can be separated into three groups as follows, A first page 
reiterates the input, thereby providing the program user a permanent record of 
the problem actually submitted to the computer. The second page, or group of 
pages, provides the Cartesian and spherical magnitudes and phases of the com- 
puted field quantities. The last page provides a summary of magnitudes of the 
Cartesian and spherical fields normalized about the maximum value which is set 
to zero decibels, arbitrarily. 


COMPUTER REQUIREMENTS 

The program is written entirely in FORTRAN for the IBM 360 Mod 9i/95 
systems at the Goddard Space Flight Center. It requires approximately 600 K 
bytes of core. The running time of a job depends on the type of problem, the 
number of field points, reflector size, sampling LI, operating frequency, and 
number of feeds. 

For a typical problem with one field point, 30-foot reflector, LI = 1.0, 
frequency = 8.0 GHz, and a single feed, the computer time is approximately 30 
seconds. 


CRITICAL PROGRAM AREAS 

This computer program has been correlated with well-known approximation 
functions, antenna slide rules, and measurements over a period of more than two 
years. Computed results agree favorably with those cases which allowed com- 
parison in terms of beamwidth, side-lobe level, null position, and gain. Never- 
theless, there are constraints which apply to the generation of radiation patterns 
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when using the program. Broadly stated, the Kirchhoff method fails when the size of 
the scattering object, its radius of curvature, and that of the wave, and the range of 
observation are small in terms of a wavelength. The failure is gradual. That is, when 
D, R c , and r are 100 \ results are excellent. For 10 \ the solution is still use- 
ful, and at 1 \ the solution is probably valueless. 

As indicated previously, the basis for the formulation assumes continuously 
varying functions. Increased departure from continuity will result in gradual 
degradation of the solution. The adequacy of the sampling interval should be 
verified occasionally with increasing $ . After the main beam and first sidelobe 
have been obtained the interval should be reduced to about half its original value. 


COMPUTED RESULTS FOR ATS-F 

A selected number of pattern cuts requested by the ATS-F Project will be 
presented here. The gain degradation of the paneled design-antenna, without the 
thermal distortion encountered in orbit, is given as a function of frequency in 
Figure 2. Degradation varies with frequency since a fixed physical departure of 
the paneled structure from a paraboloid becomes more significant (in terms of 
wavelengths) as frequency increases. A summary of gain degradation versus 
frequency is included with Figure 2. Input parameters and ancillary relation- 
ships are given below. 

Input parameters to IBM 360/91: 

focal length of parabolic cylinder F = 13.18180 ft. 
feed position for paneled antenna F C; = +13.2 ft. 

Sampling LI = 1.0 

diameter of reflector D = 30.0 ft. 

number of panels N* = 48.0 

exponential N = 1.31250 in 3 = cos N @ (-10 db edge taper) 
panel half-angle a = 3.75 degrees 

Relationships used in computations 

F r = F c /cos 2 a, where F R equals rib focal-length 

F = F_ sin 2 tt /N* _ i z i f or structure without unique focal-point. 
m R 2 tt/N' 11 



Table X 

Paneled- Antenna Verification 1 on Axis (?/ b 0°) 


Integration Net 

Method 2 

Sampling (LI) 

% (db) 

2 GHz 

(x,y) Cartesian 

Parabolic 

1.0 

-81.487 

(x,y) Cartesian 

Cylinder Panels 
RFS 

1.0 

-81.475 

(x,y) Cartesian 

RFS 

0.5 

-81.486 

(x,y) Cartesian 

AFD 

1,0 

-81.475 

(x,y) Cartesian 

AFD 

0.5 

-81.486 

(cr,£) Cylindrical 

Parabolic 

1.0 

-81.474 

3 

(a-, 0 Cylindrical 

Cylinder Panels 
RM 

1.0 

-81.474 

» / 

8.15 GHz 
(x,y) Cartesian 

T . 

Parabolic 

1.0 

-81.704 

(x,y) Cartesian 

Cylinder Panels 
RFS 

1.0 

-81.701 

(x,y) Cartesian 

AFD 

1.0 

-81.701 

(cr,£) Cylindrical 

Parabolic 

1.0 

-81.685 

(cr, £) Cylindrical • 

Cylinder Panels 
RM 

1.0 

i 

-81.685 


Ideal Paraboloid reference Level = *-81,463 db. 

2 The AFD method utilizes tangents at opposite sides of a panel and a fabric line to obtain the 
normal, while the RF$ method is the 3-point,, 2-tangent approach. The RM method is similar to 

the AFD method in that normals are obtained by an interpolation process. 

3 

[cr, C ) option leads to excess areas beyond panels 

The preceding results were obtained by several of the analytical/numerical 
techniques described previously. Table 1, above, indicates the results of some 
of the comparisons. 

Distortion of the 30-f't reflector in orbit as a function of the reflector atti- 
tude with respect to the sum has been computed by others for each position in the 
24 hour period orbit. Two representative cases are hours 11 and 4. Illumination 
edge taper s were set at -10 db, and linear polarization was used in all the follow- 
ing computations. 

The behavior of the paneled antenna under thermal load at Hour 11 and 
Hour 4 is now given at 2.0 GHz and 8.15 GHz. At each frequency a pattern cut 
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Figure 2. Gain degradation A G vs frequency for 30-ft ATS*F antenna 


on an ideal focal-point paraboloid and an undistorted, optimally-fed paneled 
antenna are provided to serve as a basis for comparison. Figures 3 and 4 are 
the paraboloid and paneled-reflector reference patterns at 2.0 GHz, while 
Figures 5 and 6 show results at Hour 11 and Hour 4, respectively, for that 
frequency. It can be seen that the paneled antenna is a very good approximation 
to a paraboloid at 2.0 GHz. The principal differences in computed results are 
a slight perturbation of the phase plot and increased level of the cross -polarized 
energy due to the panel effect. Cross -polarized levels for the ideal paraboloid 
were extremely low and are not plotted. 

At Hour 11 and Hour 4 the reflector performance at 2 GHz is seen to be 
affected only slightly. A part of a planar cut along the line of symmetry for the 
thermal distortion is given in Figures 5 and 6 to illustrate beam widening (4%) 
and axial shift due to reflector distortion. Further effects on phase character- 
istics and cross -polarization levels are in evidence but these are realtively 
.Ainor perturbations of the beam characteristics. 
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Figure 4. Optimallyfed paneled antenna (2.0 GHz) 
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Figure 6. Paneled antenna with HOUR-4 distortion (2.0 GHz) 



Proceeding to the higher frequency, 8.15 GHz, Figures 7 and 8 are the para- 
boloid and paneled-reflector reference patterns while Figures 9 and 10 show re- 
sults at Hour 11 and Hour 4, respectively, at this frequency. The paneled antenna 
is still a very good approximation to a paraboloid at 8.15 GHz. For the pattern 
cut shown, negligible beam broadening appears. The phase characteristic is now 
highly asymmetric, as is the cross -polarization pattern. Since the ideal para- 
boloid cross -polarization levels at 8.15 GHz were reported as about -80 db rela- 
tive to the principal polarization E 0mlx these were not plotted. Null-filling 
is seen to be more severe at 8.15 GHz than at 2.0 GHz. 

Table 2, below, indicates the results of some comparisons made among the 
different analytical/numerical techniques for the thermally distorted surfaces at 
Hour 11 and Hour 4. 


Table 2 


Thermally-Distorted Paneled-Antenna Verification on Axis ( 6 = 0°) 


Integration Net 

Method 

Sampling (LI) 

E * < db > 

AG = G - G . . . . , 

paneled distorted- 

paneled 

(db) 

2 GHz, Hour 11 





(x,y) Cartesian 

RFS 

0.5 

-81.495 

0.009 

(x,y) Cartesian 

AFD 

0.5 

-81.498 

0.013 

(x,y) Cartesian 

AFD 

0.25 

-81.502 

0.017 

(at 0 Cylindrical 

RM 

0.5 

-81.490 

0.016 

8.15 GHz, Hour 11 





(x,y) Cartesian 

RFS 

1.0 

-81.903 

0.202 

(x,y) Cartesian 

AFD 

1.0 

-81.908 

0.207 

(x,y) Cartesian 

AFD 

0.5 

-81.913 

0.212 

(cr, l) Cylindrical 

RM 

1.0 

-81.896 

! 

0.211 

2 GHz, Hour 4 

| 




(x,y) Cartesian 

RFS j 

0.5 

-81.496 

0.010 

(x,y) Cartesian 

AFD 

0.5 

-81.501 

0.015 

8.15 GHz, Hour 4 





(x,y) Cartesian 

RFS 

1.0 

-81.931 

0.230 

(x,y) Cartesian 

AFD 

1.0 

-81.939 

0.238 
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POLARIZATION 




Additional detailed comparisons were made, for the normals to the surfaces by 
different methods for example, but these are not presented here. Comparison 
of the general methods, for the special case of zero distortion, with the ele- 
mentary parabolic cylinder panel approach has served to establish the validity 
of the general methods. 


SUMMARY 

This document describes the analytical and numerical methods that were 
employed in obtaining the far -field or Fraunhofer radiation patterns for the 
ATS-F paneled reflector antenna via vector Kirchhoff theory. The emphasis 
throughout has been on using analytical methods whenever practicable, and 
numerical methods for the description of distorted surfaces and subsequent 
integration. A gradual transition was made from paraboloid to a paneled an- 
tenna, composed of portions of a parabolic cylinder, and finally to a distorted 
paneled antenna. By this means the analysis is justified, and numerical verifi- 
cation can be obtained at various levels of the development. For example, the 
components of the normals, the increments of area, and finally the computed 
field values can be compared in the transition. 

Although this task was undertaken to obtain solutions for a specific ATS-F 
geometry (D = 30 ft., N = 48 panels, f = 2.0, 8.15 GHz, etc.) the program is gen- 
eral in nature and can accept parameter changes without any modification. The 
assumption of straight fabric lines is not essential, and the methods for dealing 
with doubly curved surfaces are available. At this time, the more general option 
has not been implemented since it is not known whether adequate coordinate data 
or empirical curves between supporting ribs will be used to describe the surface 
to the resolution required by the sampling theorem. Feed geometry can also be 
changed to include monopulse configurations, line feeds, etc. without any modifi- 
cation of the basic program. 
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APPENDIX A 


THE INTERSECTION OF A PARABOLIC 
CYLINDER AND A PLANE 

Parametric equations for a paraboloid (F): 

. _ <7* 

x = o sin y = - cr cos £. z = — ♦ z. 

4 F 1 

Parametric equations for a parabolic cylinder F c ): when 

£ = 90°. z j = 0: 

. ^ 2 

x = (7, y = any value, z ■ £ — 

4F f 

2 

,\ z = — for a parabolic cylinder (for all values of y). 

< F e 
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equation of the intersecting plane (FI): 


r * v, s a rr x v +yv + zv 

1 x r J y x z 


r « generic vector on n 

Vj = unit vector orthogonal ton - v i + v J + v k and (v 2 + v 2 + v 2 ) 1/2 ~ 1 

• y ^ x y z 

a - minimum distance toll , a scalar 

A simultaneous solution of the equations for parabolic cylinder and plane 
yields the equation of their intersection. If v z a 0, then n contains the displaced 
z-axis through vertex point V as shown above. The intersection assumed here 
is then not of the most general type. Substituting x from the equation for the plane 
into the equation for the parabolic cylinder gives: 


z 



2 

for 


the intersection (fl). 


The singularity at v x = 0 poses no problem since it represents the trivial case 
of the intersection between the two surfaces. 

Now the displacement of V t* <m O is 


P e - 1 Xj + J yj + k z t = j y t = y t , 


above, and 


Pr • v , = x, V V + y, V v + Z, V = y, v , = a. 


The arbitrary separation y x between 0 and V can be eliminated by letting y = y x 
+ Ay so that 
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z = _L 

"Vj v y - (y, + A y) v y - z v f “ 

2 1 

Ayv y " 

4F c 

V X J 

'4F C 

L v x J 


for n. 


But tan a = v x /v , where a is the angle n makes with the displaced xz-plane. 
Also, Ay = x' sin a where x' is the distance from a point on n to the displaced 
z-axis. Then z = x’2/4(F /cos 2 a ) for the intersection, which represents a 
parabolic arc lying in plane n . The focal-length of this parabolic arc is F R = 
F c /cos 2 a. For a reflector composed of indexed portions of a parabolic 
cylinder F c , a becomes the panel half-angle. 
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APPENDIX B 

ILLUMINATION EDGE-TAPER FOR FOCAL- POINT -FED PARABOLOIDS 


The parametric -equations for a paraboloid with vector at - -F are 


.q*2 

X a or sin £, y ™ - O- cos £, z = + z. - - F 

4 F 1 4 F 


Since y . is symmetric about z, take £ -ti / 2 without loss of generality. The 
field at vertex V is E t « 3 0 /F. The field at some edge (o-) is 
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d b = 20 


+ 1 + 20 N log, 


f A F 2 


(space taper) 


(feed taper) 


In the above, db > 0, and the equation can be used to determine N for a specified 
edge taper or, conversely, to determine the edge taper if directivity is known. 


v 

i 
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APPENDIX C 


POWER INTERCEPTED BY REFLECTOR USING 3 cos N © 

Since - v < N <+«j , the effect of N on the energy incident upon a reflector is 
of interest. For 0 < N < + », the maximum value of the feed function is always 
unity (bounded above) and occurs at © = 0. Interest is restricted to the half- 
space 0 < H < 77 /2 in the present discussion 1 to avoid the generation of complex 
values for R n . For - x < N < 0, the minimum value of the feed function is al- 
ways unity (botmded below) and occurs at 0 =0. For N = 0, an isotrope (in 
amplitude) is obtained. See Figure 13. 

The use of 3 0 » cos N @ in a computer program therefore directs vary- 
ing amounts of energy at the reflector, depending on the value of N chosen. For 
some studies this is undesirable, and a constant energy to the reflector is 
preferred. A simple integration over the appropriate solid-angle allows compu- 
tation of coefficient(s) S, such that 3 0 - cos N 0 will direct equal energy at the 
antenna surface (single focal-point feed, here). In general, total energy $ r is 
proportional to 


r 0 , 

i - 

J (h> _ 


2 7T 


S 2 cos 2N 0 sin 0 d 0 d 0 = - 


= 0 *'$= 0 


2 77 S 2 
2 N + 1 


cos 


2N + 1 


0 


0, 


2 77 S* , , on + 1 

_ (1 - COS^ Ntl 

- 2N + 1 



0 < 0 2 < 77/2. 


Suppose that S = 1 (reference) and N = m for a first case and N = n for a second 
case. Equating energies leads to 

s 2 n , + 1 (1 ~ cos 2m+1 0 2 ) 1/2 

2m+l (l_cos 2n + 1 0 2 ) 


the required coefficient for the second case. 


1 The FORTRAN program actually forms 3 0= | cos © | N , which effectively avoids complex 
values for 3 0. If several feeds are employed, initial phase is introduced independently for 
the array via in 2 S . cos N * 0 . e * i^i, The energy over a solid angle can still be obtained, 
but not as simply as above. 
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Figure 13 . Prime-feed function 3 







